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Measurements of the differential cross-section and the transverse single-spin asymmetry, An, vs. 
if for 7T° and rj mesons are reported for 0.4 < xf < 0.75 at an average pseudorapidity of 3.68. A 
data sample of approximately 6.3 pb _1 was analyzed, which was recorded during p 4 + p collisions 
at y^s = 200 GeV by the STAR experiment at RHIC. The average transverse beam polarization 



3 



was 56%. The cross-section for 7T° is consistent with a perturbative QCD prediction, and the rj/n 
cross-section ratio agrees with previous mid-rapidity measurements. For 0.55 < xf < 0.75, An for 
rj (0.210 ± 0.056) is 2.2 standard deviations larger than A N for tt° (0.081 ± 0.016). 



A well known prediction of collinearly factorized per- 
turbative Quantum Chromodynamics (pQCD) is that the 
cross-section for forward meson production in proton- 
proton collisions should have negligible dependence on 
the transverse polarization of the incident proton 
This early prediction was contradicted by measurements 
of sizable pion transverse single-spin asymmetries 
(An), defined for a forward moving polarized beam scat- 
tering to the left and with a vertical spin quantization 
axis as 

A N = t , i - (1) 

In order to explain the large asymmetries, several ex- 
tensions of the pQCD collinear framework have been pro- 
posed. These approaches take into account the possi- 
ble spin-dependent transverse components of parton mo- 
mentum (Sivers effect 0]), the possible spin-dependent 
fragmentation of a scattered polarized parton (Collins 
effect Q), or higher- twist effects where transverse mo- 
menta related to the previous approaches are included 
in the hard scattering term of a collinear calculation 0- 
11] . A wide range of high energy polarized experiments, 
both nucleon-nucleon [12- 14] and lepton-nucleon 

EMI, 

have been performed to characterize the kinematic and 
process dependences of the asymmetries, and in the case 
of the latter, to directly test these approaches. 

For more than 20 years, we have known that the trans- 
verse asymmetries in forward pion production depend 
critically on the isospin projection (I3) of the produced 
mesons relative to that of the parent hadron. In proton 
scattering experiments, the asymmetry for the 7r~ me- 
son, which contains a down quark and an anti-up quark, 
has the opposite sign relative to the asymmetries for the 
7T+ and 7T° mesons, produced from the predominant up 
quarks. 

In this paper, we report for the first time at y/s = 
200 GeV the transverse single-spin asymmetry for the rj 
meson, another member of the pseudo-scalar octet that 
has the same isospin projection as the it (I3 = 0). We 
note that the FNAL-E704 collaboration previously found 
a large An for the 77 for Feynman-x (longitudinal momen- 
tum of the observed particle divided by the beam energy) 
> 0.4 at y/s = 19.4 GeV jl^]. In addition, we report the 
differential cross-section for rj production in the region 
where the spin asymmetry is measured. 

Leading-twist collinear pQCD has been successful in 
describing a wide range of unpolarized cross-section mea- 
surements at the Relativistic Heavy Ion Collider (RHIC), 



from 



at forward rapidity [20] l2j, to tt° [22J-[2I 
rj lH-izt tt ± /K ± '28\, and jets [I| at mid-rapidity. 
Such agreement is considered a strong indicator that the 
given process can be interpreted within the framework of 
pQCD. Therefore, the comparison between the unpolar- 
ized cross-section and the leading-twist collinear pQCD 
prediction becomes the basis on which to apply the afore- 
mentioned theoretical extensions to the associated trans- 
verse spin effects. 

For forward 7r° production, recent STAR measure- 
ments of the cross-section are consistent with next-to- 
leading-order (NLO) pQCD calculations in the same re- 
gion where a large transverse spin asymmetry is found 
20, 21]. However, these results do not cover the large 
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Feynman-x (xf) region where the acceptance for the rj 
decaying into two photons becomes large. In this paper, 
we have extended the analysis of the ir° cross-section and 
An to xp of 0.75, where its spin asymmetry and cross- 
section can be directly compared to those of the 77 mesons. 

The data were taken with the STAR Forward Pion De- 
tector (FPD). The FPD is a modular lead glass calorime- 
ter located at forward rapidity in the STAR interac- 
tion region at RHIC at Brookhaven National Labora- 
tory. Two modules were placed on either side of the 
beam line, covering the pseudorapidity region from ap- 
proximately 3.3 to 4.0. Each module contained 49 cells 
(glass blocks approximately 18 radiation lengths deep), 
forming a 7 x 7 square array. The data were collected 
in 2006 with transversely polarized proton beams and 
an integrated luminosity of ~ 6.3 pb _1 . The average 
polarization was (56.0 ±2.6)% for the beam facing the 
FPD. As An is a single spin observable, the spin state 
of the second beam (with (55.0 ±2.6)% polarization) was 
integrated over for the An at positive xf- At negative 
x f , the spin state of the first beam was integrated over. 
Events were recorded only when the total ADC count in 
either of the two modules was greater than a fixed thresh- 
old, which was nominally equivalent to 30 GeV. Photons 
reconstructed within a quarter of a cell from the detector 
edge were discarded. Only those events with exactly two 
reconstructed photons were analyzed, with the resulting 
loss of yield corrected for the cross-section measurement. 
The STAR Beam Beam Counter (BBC) on the away side 
was used to reject the single-beam background. The near 
side BBC was not required to produce a signal, as most of 
the analyzed events already had more than half the beam 
energy deposited in the FPD. The efficiency for the away 
side BBC condition was estimated to be (93 ±4)% for all 
non-singly diffractive events based on previous analyses 

The xf coverage of the previous analysis of Ajv(7r°) 
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FIG. 1. (color online) The distribution of 0\l O9 , as denned in 
the text, for Eduster > 65 GeV for data and simulation, in 
units of FPD transverse cell size (3.81 cm). For comparison, 
the one and two photon cluster peaks from simulation were 
independently normalized, and uniformly shifted by +0.01 
transverse cell size (1 bin) to account for the small difference 
in the average size of clusters between simulation and data. 
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J12j, which included this data set, was limited by the 
difficulty in separating 7r° clusters from single photons 
with xf > 0.55. At this point, the typical separation of 
two 7T° decay photons at the surface of the FPD becomes 
similar to the Moliere radius of the lead glass (3.32 cm) 
and transverse cell size (3.81 cm). On the other hand, 
the rj meson acceptance lies mostly above an if of 0.5 
due to the larger separation of its decay photons. 

With the current analysis, the tt° - 7 separation has 
been greatly improved by analyzing a Log, defined as 



Oho 



'Ei Log[(Ej + E )/GeV] ■ (x - Xl f 
Log[(E i + E )/GeV] 



(2) 



where Ei and Xi are the energy and the location of the 
i th channel, and Eq = 0.5 GeV. The i th term in the sum 
is skipped if Log(Ei + Eq) < 0. It provides a significant 
sensitivity to the topological differences between single 
and double photon clusters at high energies, as evidenced 
by the clear separation between the one and two photon 
peaks in Fig. [1] Also shown are the results from PYTHIA 
and GEANT simulations, which closely reproduce the <7Log 
distributions for both types of clusters up to a small offset 
(~ 1% of the transverse cell size). As a result, the xf 
coverage for 7r°'s was extended from 0.55 to 0.75. 

In addition, the GEANT simulation of the electromag- 
netic shower in the FPD is now based on the tracking of 
optical photons produced by the Cherenkov effect. Com- 
pared to the previous method based on charged particle 
energy loss, the new simulation produces a better agree- 
ment with the data on shower shape, energy resolution, 
and the observed shift in gain as a function of photon 
energy. Combined with a more advanced parameteriza- 
tion of the shower shape including the effects of incident 



FIG. 2. (color online) (a) Di-photon invariant mass, M 77 , 
distributions in data and simulation for _E 77 > 45 GeV, with 
the "center cut" as defined in Eq. (0 . The simulation results 
were normalized to have the same number of events as the 
data in the vr mass region (0.08 < M 77 < 0.19 GeV/c 2 ). The 
symbol (rf) indicates the average pseudorapidity of the photon 
pair, (b) Same as (a), but plotted using an expanded linear 
scale to illustrate the r\ mass region. For the dashed line, the 
r\ signal was removed from the simulation at the pythia level, 
(c) 

bars are statistical uncertainties only. 



angle, it allows for a higher precision calibration needed 
for the cross-section measurement. 

The top two panels of Fig.[2]show data-simulation com- 
parisons of the di-photon invariant mass spectra. The 
"center cut" , so named because it covers roughly the cen- 
tral region of the FPD acceptance, is imposed on all event 
samples in order to enhance the r\ meson acceptance rel- 
ative to the background. It is defined as 



(f? 77 - 3.65) 2 + tan 2 



"77. 



< 0.15, 



(3) 



where ry 77 is the pseudorapidity of the di-photon center 
of mass relative to the polarized beam, and c/> 77 is its az- 
imuthal angle. The distributions of ir° and rj events in 
the FPD, and the subset of each that pass the center cut 
are shown in Fig. [3J A full simulation based on pythia 
6.222 and GEANT 3 was compared to the data. The reflec- 
tivity and absorption properties of the aluminized mylar 
wrapped glass blocks were varied to minimize the dis- 
crepancies between the photon shower shape in the sim- 
ulation and that measured in the data. While detailed 
knowledge of the glass-mylar interface remains a limit- 
ing factor in the precise modeling of the shower devel- 
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STAR p + p -> 7i°, r| + X at sfs = 200 GeV 
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tan(*„) 



FIG. 3. Pseudorapidity vs. tangent of the azimuthal angle 
of the di-photon center of mass, for i? 77 > 50 GeV. LEFT: 
0.08 < M T7 < 0.19 GeV/c 2 , RIGHT: 0.45 < A/ 77 < 0.65 
GeV/c 2 . The filled boxes indicate events that pass the center 
cut (Eq. ©). 



opment, the agreement in the widths of mass peaks be- 
tween the simulation and data has been improved signif- 
icantly over previous analyses 12, 2(| 21 1. Furthermore, 
the data-simulation agreement in the continuum region 
between the 7r° and r] peaks is very good, allowing for a 
simulation-based background estimation for the rj signal. 
Corrections for the remaining data-simulation discrepan- 
cies in mass resolution were applied to the cross-section 
measurements. The ?y to 7r° cross-section ratio in the 
simulation has been set at 0.45 to be consistent with the 
data. The bottom panel of Fig. [5] shows the invariant 
mass dependence oi An, which exhibits a suppression in 
the continuum region. Within the large statistical uncer- 
tainty, the asymmetry for this region does not show a sig- 
nificant xp dependence. In the simulation, this mass re- 
gion is dominated by approximately equal contributions 
from a pair of photons from two different ir° decays, and 
a charged hadron combined with a photon. 

The energy resolution of the FPD is estimated to be 
about 7 to 8% of the total energy based on the com- 
parison of invariant mass and di-photon separation dis- 
tributions between data and Cherenkov shower simula- 
tion. Coupled to the rapidly falling cross-section in en- 
ergy, more than half of events in any measured energy 
bin originate from lower true energy bins. For the cross- 
section measurements, we unfolded the energ y sm earing 
by applying the Bayesian iterative method [3l| to the 
smearing matrices obtained from the simulation. The 
unfolding procedure combines the statistical and system- 
atic uncertainties from the original data points. 

The upper panel of Fig. 2] shows the differential cross- 
sections for 7T° and rj. The center cut (Eq. ([3])) was 
imposed on both mesons. Full pythia + GEANT simu- 
lations were used to obtain the detector efficiency cor- 
rections including the rj — » 27 branching ratio. Also 
shown are the previously published STAR results for the 
7T° cross-section in similar kinematic regions. The error 
band corresponds to the NLO pQCD theory prediction 
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FIG. 4. (color online) Differential production cross-sections 
for 7T° and rj at average pseudorapidity of 3.68. Also shown are 
the previously published STAR results for similar kinematics 
[2lj and a NLO pQCD calculation of the n° cross-section [3^ ]. 
The error band represents the uncertainty in the calculation 
due to scale variations. The rj to tv° cross-section ratio is 
shown in the bottom panel. The error bars indicate the total 
statistical and systematic uncertainties. 



for the 7T° cross-section [32fl. based on the CTEQ6M5 par- 
ton distribution function [33| and the DSS fragmentation 
function 34|. The uncertainty for the theory prediction 
was obtained by increasing the factorization and renor- 
malization scales from fj, = pr to fi = 2pr- We note 
that the DSS fragmentation function includes in the fit 
the previously published STAR results at pseudorapid- 
ity of 3.3 and 3.8 [2(j, along with other RHIC results. 
The error bars include both statistical and systematic 
uncertainties. The major sources of systematic uncer- 
tainties are the absolute energy calibration uncertainty 
of 3%, which dominates the tt cross-section, and the un- 
certainty from the unfolding process, which dominates 
the rj cross-section at high energies. The normalization 
uncertainty was estimated at 12.5%, including the un- 
certainty of the BBC coincidence cross-section of 7.6% 

The lower panel of Fig. 0] shows the rj to ir° cross- 
section ratio, which is found to be around 50%. The 
error bars include both statistical and systematic uncer- 
tainties. The latter is dominated by the 1.5% relative 
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FIG. 5. (color online) An vs. xf at average pseudorapidity 
of 3.68 for 7T° and rj. Also shown are the previously published 
results for tv° at lower if, derived from the same data set 
but without the center cut [l^] . The error bars are statistical 
uncertainties only. The error boxes indicate the systematic 
uncertainties. 



energy scale uncertainty, caused by the acceptance dif- 
ferences for 7T° and rj decay photons, and the localized 
variations in cell to cell calibration. The absolute cali- 
bration is common to both mesons, and largely cancels 
for the ratio. 

In pQCD, large- x f production of both 7r° and 77 arises 
from hard-scattered partons fragmenting into mesons 
with large momentum fraction z (ratio of hadron mo- 
mentum to the momentum of its parent parton). The 
fragmentation process generally does not depend on the 
details of the hard scattering, and a single set of pion frag- 
mentation functions explains a wide variety of RHIC data 
H [H, HU. While there are currently no NLO pQCD 



predictions for the forward r\ production cross-section, 
we note that our measurement of the ir° cross-section is 
consistent with the NLO prediction, and the 77/71"° cross- 
section ratio is consistent with the recent NLO pQCD 
extraction of the 77 fragmentation function [3fj| . 

Figure O shows the A n, calculated using the "cross ra- 
tio" formula [13, 37] , as a function of xp for n° and 77 after 
correcting for the underlying background. Also shown is 
the previous STAR result [12[ for An(tt°) at lower Xf, 
which utilized the same data set as the current analysis 
but without the center cut. The two ir° results are con- 
sistent within their correlated errors. The background 
correction, which only significantly affects the rj asym- 
metry at medium energy, is obtained from a simulation 
sample corrected for the 77 yield and mass resolution, and 
the assumed background Am of 0.005 ± 0.0f6 extracted 
from Fig. [SJc) . The error bars indicate statistical uncer- 



tainties only, while the error boxes indicate the system- 
atic uncertainties. The main source of the systematic 
uncertainty is the background correction; the polariza- 
tion uncertainty is negligible in comparison. The An for 
negative xf is consistent with zero for both mesons. 

The comparison of An for 7r° and 77 mesons is of par- 
ticular interest given their similar up and down quark 
content, with wave- functions of both mesons containing 
uu and dd pairs. The 77 differs from the 7r° mainly in 
that it is in an isospin singlet state, and that it contains 
ss in the wave function. The latter results in 77 being 
significantly more massive than the 7r°. 

In conclusion, STAR has measured the x f dependences 
of the cross-section and transverse single-spin asymme- 
tries for 7T° and 77 mesons produced at an average pseudo- 
rapidity of 3.68 in yfs = 200 GeV polarized proton colli- 
sions. For 0.55 < x F < 0.75, A N (rj) (0.210± 0.056) is 2.2 
standard deviations larger than An(tt°) (0.08f ± 0.0f6). 
In this kinematic region, both the ir° cross-section and 
the 77/77° cross-section ratio are consistent with NLO 
pQCD expectations. This suggests that the measured 
77 asymmetry can be understood within the framework 
of pQCD. While several calculations exist for pion and 
kaon asymmetries ES mum, the first pQCD calcula- 
tion of An for the 77 meson was performed only recently 
[42| . This model generates an 77 asymmetry that is sub- 
stantially larger than that for the ir° via a sizable initial- 
state twist-3 effect for strange quarks. The calculated 77 
asymmetry rises to about 12% a.t xf of 0.4, well above 
our measured asymmetry, but then agrees quantitatively 
with the data for xf > 0.5. It is yet unknown if a sim- 
ilar difference can arise from the fragmentation process 
via the Collins effect. A higher statistics measurement 
of the An for the 77 meson in this kinematic region is 
necessary to make a precise comparison to that for the 
7r°. Understanding the exact nature of these asymmetries 
can be further aided by complementary measurements of 
An for final states that lack Collins contributions, such 
as jets and prompt photons. 
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